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1. Introduction 

In a previous series of papers P-|2], we have presented the exact diagonalization of the 
hermitian Schrodinger operator corresponding to a periodic potential composed of 
atoms modelled by delta functions with different couplings for arbitrary A^. This basic 
structure can repeat itself an infinite number of times giving rise to a periodic structure 
representing a A^-species one dimensional infinite chain of atoms. To our surprise this 
model is exactly solvable and a far from straightforward calculation leads to an exact 
band condition. Due to the factorizable form of the solution one has not only the 
advantage of closed form expressions but one can also perform computer calculations 
with an exceedingly degree of accuracy due to the exact nature of the solution itself. 
The main physical motivation was indeed to modelate the band structure of a one 
dimensional quantum wire. The model can be extended with a minimum amount of 
effort to the non-hermitian but PT-symmetric [3| quantum hamiltonian and the band 
condition becomes real although the potential is obviously a complex one. The imaginary 
parts of the couplings play a central role in defining the band structure which can be 
modeled at will by varying this parameter. The purpose of this paper is to present 
this new solution in a detailed fashion. In a previous paper a preliminary analysis of 
the existence and properties of this exact band model solution were given j^. Here we 
want to present a careful study of the solution in a more rigurous manner paying due 
attention to the exciting properties that this rich structure seems to exhibit. 

The idea of performing band structure calculations by using a non-hermitian 
Hamiltonian is extremely promissing as the examples of PT-symmetric quantum 
hamiltonians so far existing in the literature rely more in aspects concerning the discrete 
spectrum and bound states ([Sj, [Oj, [Z| and |H]) and are also far from an actual physical 
application. It is not the aim of the present paper to provide a full discussion of 
the manifold aspects of the PT-symmetry in quantum mechanics. We address the 
interest reader to the original reference |3j , some recent interesting theoretical criticism 
jHj and also to some work done recently in one dimensional models which may help 
to understand better the role of PT-symmetry in the framework of one dimensional 
systems. In all these papers the imaginary part of the coupling plays a substantial role 
in defining the spectrum of bound states. We shall show that in the case of the band 
structure hereby presented this imaginary parts are of primary importance in defining 
the transport properties of the one dimensional quantum chain. 

The plan of the paper is the following. In Section |21 we shall be discussing the 
main features of the model just as a simple model for a quantum wire with only real 
couplings. Furthermore the presence of "PT-symmetry and the role of the complex 
couplings will be introduced. The band spectrum remains analytical and real in spite of 
the complex non-hermitian nature of the potential. In Section El a throughout analysis 
of the band condition is made in order to clarify the effect of the imaginary parts of 
the couplings on the band spectrum. This analysis can be expressed analytically for 
some cases where N is not too large. Otherwise the equations appear to be intractable. 
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Even in this case numerical calculations can always be performed with a large degree 
of accuracy. In section 0] we make some comments on the band spectrum and observe 
the form of the PT-symmetric electronic states. Several pathologies plaguing the band 
spectrum for the complex coupling case are observed. It can be avoided working on a 
small range of energies for certain values of the imaginary couplings. The localization of 
the wave function over the primitive cell is also found to be altered manipulating only 
the imaginary parameters. We close with a section of Conclusions. 



2. The model and its band spectrum 

Let us begin with a brief remainder of the solution presented in P] corresponding to a 
infinite periodic potential composed of a basic structure made out of a finite number 
N of equally spaced deltas with different real N couplings, repeating itself an infinite 
number of times. If the spacing is a and 

aj = (2.1) 

is the length associated to each coupling (ej), the band structure can be written as 

cos{NQa) = B{e]ai,a2, . ■ . ,aj\f) (2.2) 
where Q is an arbitrary real number which varies between and —-^ and 

B{6; a,,..., a^) = 2^-^ ^ h,...{N)...hk - 2^^' Yl ^-(^ " 2).../i, 

p p 

+ 2^"^ J2h^...{N- 4).../ifc - . . . (-l)f -^2 J2 K...{2)...h, 
p p 

+ (-l)f (2.3) 

B{e- ai, . . . , aN) = 2^"^ J] h,...{N)...hu - 2^-^ ^ h,...{N - 2)...hu 

p p 

+ 2^-^ ^-(^ -A)...h,-... (-1)^2^ h^..m...h, 



p 

N-l 



+ (-!) — (/ii + /i2 + ... + /ijv) (2.4) 

for even and odd repectively. 

The symbol J2p hi...{M)...hk means a sum over all products of M different hi 's with 
the following rule for each product: the indices must follow an increasing order and to 
an odd index must always follow an even index and reciprocally. 

The functions hj have the universal form: 

, / \ / \ f \ sin(£:) , ^, 

hj{e) = cos{e) + j^— J (2.5) 

and the independent variable is a function of the energy (i.e. e = ka). In order to see 
that this condition looks much simpler than one might think at the beginning of the 
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calculation let us list below, for the benefit of the reader, the first three conditions for 
the cases = 2, 3 and 4. 

cos(2Qa) = 2hih2 - 1 (2.6) 
cos(3Qa) = 4/ii/i2/i3 - + /la + /^a) (2.7) 
cos(4Qa) = 8/11/12/^3/^4 - 2(/ii/i2 + hih^ + /la/^s + h^h^) + 1 (2.8) 

It does not require too much time to write down the band conditions for fairly large N, 
but more important is the fact that the exact formulae ()2.3|) and ()2.4|) are in itself quite 
easy to program for sequential calculations. Using the band condition we have been able 
to generate not only a full band spectrum but also to calculate the density of states for 
periodic one dimensional lattices with a high degree of accuracy (See references llj and 

The next step is to include PT-symmetry in this model. It is not hard to see that 
()2.2|1 is still real if we include the following changes: 

• Promote the couplings from real to complex, i.e. — >■ rj + isj 

• Order the potential in a "PT-invariant form, that allow us to choose a PT-symmetric 
primitive cell. This leads to the following identifications: 

hN = hi 
hN-i = hi 



/liv , 1 = h*N even 
h N+1 = h*N+i odd 

It is easy to check that equations (j2.3j) and (j2.4p remain real under these identifications 
which make obviously the periodic potential complexhut PT-invariant. There has been 
an earlier attempt to generate real band condition from a complex but PT-invariant 
potential ^1]. However the results concerning the appearance and disappearance of 
forbidden and allowed bands was inconclusive. In our case this effect is clear and will 
be discussed at length below. 

As is well known for years a hermitian periodic potential cannot alter its band 
spectrum just by fine tunning the couplings. The bands can indeed be made wider 
or narrower but its number and quality (forbidden or allowed) remains unchanged. 
The theorems supporting these statements are all based upon the intuitive idea that a 
hermitian operator cannot change its spectrum that is basically given by the eigenvalues 
and eigenfunctions of the states at the edges of the bands. The mathematics can be 
hard but the physical idea was indeed whether this behaviour would be mantained if a 
PT-invariant potential is used. For this purpose they use various analytical potentials 
carefully shifted to be PT-invariant. They can prove that the band condition is real 
but in order to analyze the band structure they have to assert with a very high degree 
of precision whether a given curve is above (below) +1 (—1) in a similar manner as we 
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have to ascertain ourselves that the expressions (|2.6|) - (j2.8|) (and in general (j2.3p and 
(12 .41) ) exceeds +1 or goes below —1. In the case of [llj this appears a very hard task 
indeed as the authors do not have to their disposal an analytical band condition, so they 
must carry out various kinds of approximations. The authors conclude that " despite this 
impressive precision, [...its equations] (16) and (17) cannot he used directly to answer 
the crucial question of whether there are hand gaps hecause these approximations to the 
discriminant [... hand condition] never cross the values [...normalized to] ±1" 

But we do have such an exact band condition and in spite of the apparent formidable 
aspect of the expressions ()2.3|) - (j2.4j) we can perform various kinds of exact and numerical 
analysis in order to check the dependence of the band width and the band number on the 
exceptional parameter that arises in our model: the imaginary part of the complex 
couplings. This will be the subject of next Section. 



3. The band condition for couplings with non-vanishing imaginary part 

In order to understand the effect of the imaginary part of the couplings on the band 
spectrum one must analize in detail the behaviour of the band condition. Let us 
begin with the simplest PT-symmetric chain including 2 deltas in the primitive cell 
characterized by the parameters (a/ai) = ri +isi and (a/ai)* = ri — isi. This case has 
been also carefully studied considering different distances between the deltas in reference 
[T^ . The band structure is entirely determined by the function B = 2 hih\ — 1, which 
can written as 

e . e(., = 0) + (3.1) 

The term due to the imaginary part is always positive and it has the effect of lifting up 
the band condition for all values of ri. The position and width of the allowed (forbidden) 
bands come from the intersections at B = 1,-1 which depend strongly on the position 
of the maxima and minima of the band condition. The changes on these limiting points 
can be analytically described in this case. 

From the equation dB/de = one can find the value of the oscillatory part of the 
band condition (i.e. the trigonometric functions) as a function of e and substituting 
into B with a properly choice of signs will result in an analytical form for the curves 
crossing the extremal points of the band condition which reads: 

C± = 2 - 1 (3.2) 

where C+(C_) lays on the maxima (minima) and the functions f± are defined as follows: 

/± = [%/2ri(ri - is,) ± (F± - 2rle'y] (3.3) 

F± = rl + {rl + slY + 2(ri + - sl)e^ + ± 

(ri -rl-sl + e^) [{r, + rj + slf + 2{n + rf - sl)e^ + e^]K (3.4) 




Figure 1. Band condition for iV = 2 when s ^ r and s ^ r. The coloured hnes 
corresponds to C±. 



Hence, the behaviour of the band condition can be studied through the evohition of C±. 
Taking the hmit si — >■ we obtain, 

\imC±={ ri(2 + ri) +£2 + r| (3.5) 



-1. 



The real part of the couphngs controls the amplitude of the oscillations (and therefore 
the distance between C+ and C_) but the minima will always lean on i3 = — 1 unless 
the imaginary part is non zero. Thus as we increase si the minima will rise, as can be 
seen in figure HJ^a), and they never touch —1 because lim^^oo C- = —1- On the other 
hand considering ri — > 0, after some algebra one obtains 

r ^ , 2sl si + e' + sl{l-2e') 

-1 H 2 ^01 £ < Si 

\imC+={ si + e^ + slC^-2e^) (3.6) 

[1 for £ > Si 

and for C_ the limit is the same as for C+ but interchanging with respect to the intervals 
in e. Thus for an imaginary part of the coupling large enough compared to the real one 
the band condition goes above +1 for energies below the value Si. Figure ^b) shows 
this last configuration. In fact the situation si ^ ri can be easily understood from the 
form of the band condition that can easily be written in this limit as: 

i3 ^ -1 + M + 2 cos\x) f 1 - 41 (3-7) 



£2 ' ^ M £2 



which shows clearly the boundaries when e < Si : l<i3<— 1 + 2s1/e'^ saturated at 
e = nn and e = {2n + l)f n G Z respectively. For e > Si the boundaries simply 
interchange among them. Due to the nature of the functions hj (e) the band condition is 
always tied up to the value +1 at every multiple of vr for all ri, si, as figure ^b) clearly 
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(a) ri = 4, r2 = 5, si = 3 



(b) n = 3, r2 = 2, si = 5 





(c) n = 3, r2 = 2, si = 25 



(d) n = 10, ra = 10, si = 25 



Figure 2. Band condition for = 3 for different values of the couplings. Each type 
of extremal point is labelled with a different colour. 



shows. Several configurations of tlie spectrum as tlie ones in figure C] can be built for 
different values of the complex coupling. 

Let us now consider a PT-symmetric chain whose primitive cell includes 3 atoms 
charaterized by the parameters (a/ai) = ri + isi, (0/02) = r2 and (a/ai)*. The band 
condition function reads B = 4 hih\h2 — {hi + hl + h2). This can be arranged as: 

As\ sin^(£:) [e cos{e) + r2 sin(£:)] 



B = Bisi = 0) + 



(3. 



The extra term arising from the imaginary part of the coupling can be either positive 
or negative depending on the value of r2 and the energy. The effect of Si is different 
compared to the previous case although there is only one imaginary parameter. The 
structure of the band condition seems to be composed by pieces each one including 
3 different types of extremal points, as figure Efa) shows, and each type apparently 
following a given pattern. The curves crossing the extremal points of the band condition 
can also be analytically calculated in this = 3 case but the expressions are quite hard 
to simplify and therefore we do not provide the explicit form of the equations. Three 
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curves are obtained, one for each type of extremal point. Initially when si = all the 
extremal points are outside or on the borders of the range [—1,1]. As we increase the 
imaginary part, the amplitudes of two of the three types of extremal points begin to 
decrease (green and red in figure |2tb)). At the same time as Si becomes larger the green 
and red curves get narrower with decreasing energy until e ~ si. Below this value the 
curves broaden trying to expel the extremal points of the band condition outside the 
target range (figure |2fc)). From the form of B it is clear that the even(odd) multiples of 
vr will remain fixed to 1(— 1). The efficiency of this expelling process depends upon the 
values of the real couplings as they control the amplitudes of the oscillations. One can 
stretch the maxima and minima (up and down respectively) increasing ri, r2, as shown 
in figure Efd). 

Several configurations can be obtained among the ones shown in figure |2l and 
one important feature must be emphasized: the band condition always changes 
"symetrically" with respect to the abscisa axis. That is: there is the same amount 
of positive and negative function whatever the values of the couplings are. In fact the 
curves laying on the extremal points show an exact refiection symmetry around the 
energy axis. This is in great contrast to the N = 2 case where the band condition can 
be pushed up above +1 for certain couplings. Thus when si ^ ri,r2 there will always 
remain a trace of the original permitted energy ranges in the form of "fiat" bands. 

The first chain for which two different imaginary couplings can be manipulated is 
the one with = 4 atoms inside the primitive cell: (a/ai) = ri + isi, (0/02) = r2 + is2, 
(0/02)* and (a/ai)*. In this case we have not calculated analytically the curves crossing 
the extremal points. Nevertheless a pictorial approach for some of the curves has been 
included in the graphics for a better understanding of the tendencies of change. The 
initial situation (si = S2 = 0) is a common one for a periodic chain (figure Efa)) and two 
different group of extremal points can be distinguished. As we increase the imaginary 
part of the couplings one group of extremal points starts to decrease and the band 
condition lifts up globally (figure Efb)). The effect of the two imaginary parts is quite 
similar. When both of them have comparable values the two groups of minima follow 
the same tendency moving upwards the band condition. However if one imaginary part 
grows much more than the other one, several minima strech down in a region of energy 
roughly included in Ssmaii < £ < Siarge (figure El^c)). Finally when Si,S2 ^ ri,r2 one 
can force the band condition to go above +1 (keeping our well known knots at e = rivr) 
in a certain energy range which depends on the values of the imaginary parts of the 
couplings (figure Efd)) approximately the same way as for the N = 2 case. Unlike the 
A^ = 3 example the band condition can be unbalanced to positive values with a properly 
choice of the parameters. 

The band condition has been also studied in detail for A^ = 5 and A^ = 6. Its 
behaviour gets really complex as the number of atoms grows. In figure ID and figure El 
some characteristic examples are shown. For N = 5 the band condition shows five types 
of extremal points that evolve "symmetrically" around the abcisa axis as we change the 
imaginary couplings. For A^ = 6 three types of maxima and minima appear and for 
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(c) n = 3, r2 = 2, si = 4,32 = 25 (d) n = 3, ra = 2, si = 20, S2 = 21 

Figure 3. Band condition for = 4 with different values of the couplings. The 
coloured dashed lines are intended to clarify the behaviour of different groups of 
extremal points. 

large enough values of the imaginary couplings the band condition is greater than +1. 

To summarize this mathematical analysis of the band condition let us now list some 
of the most salient features: 

• the band condition is composed of pieces integrated by N{N/2) extremal points 
for odd(even) A^. Each extremal point belongs to a group that evolves differently 
according to the values of the imaginary parts of the couplings 

• for odd the band condition is "symmetric" around the abcisa axis for all values 
of the couplings. Therefore, some permitted levels always remain as a part of the 
spectrum 

• for even the band condition can be expelled out of the target range for certain 
values of the couplings and so removing the allowed bands 

• for 6 = UTT the band condition is fixed to +1 or —1 depending on the parity of N 
and n for all values of the couplings 
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50 60 



(a) n = 4, r2 = 1, rs = 1, si = 20, sa = 5 (b) n = 4, ra = 2, = 3, si = 10, S2 = 20 



Figure 4. Band condition for a T-'T-symmetric chain with N = 5. 




(a) n = 5, r2 = 3, ra = 4, si = 2, S2 = 5, (b) n = 1, r2 = 2, rs = 0.5, si = 10, S2 = 11, 
S3 = 4 S3 = 12 



Figure 5. Band condition for a PT-symmetric chain with = 6. 



• the real parts of the couphngs are always proportional to the amplitude of the 
oscillations. 

Even with all the included figures it is hard to make a complete explanation of the 
behaviour of the band condition that otherwise can only be fully understood watching 
some proper animations of the function. This is the procedure we have followed to 
support the results. We encourage the interested reader to reproduce those animations 
that can easily be done with Mathematica. Some example code lines are listed in the 
appendix. 



4. The band structure and the electronic states 



Let us now make some comments about the band structure. One important feature 
that stands out from the above examples is the presence of maxima(minima) inside 
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(a) (b) (c) 

Figure 6. Band structure in the irreducible part of the first Brillouin zone for different 
chains: (a) = 3 with ri = 1, r2 ~ 2, si = 3. Different colours show the first four 
bands, (b) N — 5 with ri — 4, r2 — 2, = 1, si = S2 — (red) and Si = 0.775, 
S2 = O.f (blue), (c) = 6 with ri — 5, r2 — 2, — 3, si — S2 — S3 = (red) and 
si = 1.6f3, S2 = O.f 2, S3 = 0.3 (blue). 



the target range (—1,1). This fact involves the presence of points in the recripocal 
space where the gradient of the energy diverges that makes the physical interpretation 
of the hamiltonian not so straightforward. Probably some additional restrictions have 
to be imposed on the PT-symmetric hamiltonian. These special points also cause some 
changes on the form of the bands which are always understood as the ensembles of 
states and energies labelled with a certain index {n for the n-th band). Every point 
of the Brillouin zone must represent a physical state on every band. Therefore for a 
certain value Qa let us label the eigenenergies of the hamiltonian following an increasing 
order {ei{Qa) < e2{Qa) < . . .). These procedure leads unavoidably to the definition of 
bands whose dipersion relations are not continuos over the reciprocal space as seen in 
figure ini^a). Of course the choice of the bands indices is not unique but none of them 
will be free of these discontinuities. 

On the other hand those pathologies can be avoided if one restrics his working 
scenario to an small range of energies. In that case the imaginary parts of the couplings 
can be tuned to modify essentially the spectrum of the system removing or decreasing 
gaps to change virtually its response to transport phenomena (figures IHfb) andlH^c)). 

One can also wonder about the form of an electronic states belonging to such a 
characteristic spectrum. It is not hard to calculate analytically the wave function inside 
the primitive cell as a function of the position and the energy with the help of the 
computer for low A^. To our surprise we have found that the imaginary parts of the 
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Figure 7. Electronic wave function normalized to 1 inside the primitive cell of a VT- 
symmetric chain for TV = 3 with ri = 5, r2 = 3. si — O(blue) and si = 19(red) at 
e = 5.2. 

1 

0.8 
0.6 
0.4 
0.2 

12 3 4 

x/a 

Figure 8. Electronic wave function normalized to 1 inside the primitive cell of a VT- 
symmetric chain for = 4 with ri = 3, r2 = 4, si = 5, S2 = 8.1 for different energies: 
e = 2.95(red) and e = 3.13(blue). 

couplings behave as control parameters of the localization of the electrons inside the 
primitive cell. Thus one could tune these parameters to decrease the probability of 
presence in several sectors almost to zero or distribute it more homogeneously over 
the primitive cell, as shown in figure [7| Also for non- vanishing imaginary parts of the 
couplings the depence of the state on the energy seems quite strong as a small variation 
of this energy can involve an important change in the shape of the state (figure |H1). 

5. Conclusions 

In this paper the band structure of a one dimensional A^-difFerent coupling delta periodic 
potential is also assumed to be PT-symmetric and hence the delta couplings can be made 
complex. In spite of the non-hermitian nature of such a Hamiltonian, the band spectrum 
condition remains real. The imaginary parts of the couplings play an essential role in the 
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variety of shapes and forms that the band spectrum and the wavefunctions may acquire 
in this new framework. Ahhough the band condition is analytical for arbitrary A^, we 
have studied carefully just some cases with a low number of couplings. A throughout 
study of the N = 2 case is given with analytic expressions for the limiting points of 
the bands. The cases A^ = 3, 4, 5 and 6 are solved graphically and conclusive results 
for the band spectrum are also offered. In all these cases, the control of the imaginary 
parts plays an essential role in the way that the band spectrum gets distorted and 
reshaped. As we have been able to discuss the evolution of this band spectrum, we are 
now planning to apply these results to more physical situations in which localization and 
de-localization may be tuned on and off by means of a judicious control of the imaginary 
parts. We have been able to show how wavefunction de-localization takes place actually 
for particular cases as the ones discussed in Section ^ However a manifold of interest 
questions remains to be answered such as: Could a random (non-periodical) calculation 
with this sort of potentials be carried out in the case of PT-symmetry ?, Does the de- 
localization property holds for more general cases?. These and related research subjects 
are now being investigated and will be reported elsewhere. 
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Appendix 

Here are some example lines for animating the band condition. First load the package 
"Animation" and define the band conditions for N = 2 and A^ = 3. 

<< Graphics 'Animation' 

h [r_ , s_] =Cos [x] + (r+I*s) *Sin [x] /x ; 

B2 [rl_ , sl_] =2*h [rl , si] *h [rl , -si] -1 ; 

B3 [rl_ , r2_ , sl_] =4*h [r 1 , si] *h [r2 , 0] *h [r 1 , -si] 

-(h[rl,sl]+h[r2,0]+h[rl,-sl]) ; 

To see the evolution of the band condition for A^ = 3 versus the imaginary part si for 
ri = 2 and r2 = 3 one would write: 

Animate [Plot [Evaluate [{B3 [2 , 3 , t] , 1 , -1}] , {x , , 30} , 

PlotRange->{-3 , 3}] , {t , , 30}] 

and finally one must group and compress the cells with all the generated plots and by 
double clicking on the image the animation will start. 
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